Improved Asteroid Astrometry and Photometry 

with Trail Fitting 

Peter Veres 1 , Robert Jedicke 1 , Larry Denneau 1 and Richard Wainscoat 1 

(n: 

Institute for Astronomy, University of Hawaii at Manoa, Honolulu, 96822 HI 

<N 

P-i" and 
00 

t""*- ! Matthew J. Holman 2 

(N 

Harvard-Smithsonian Center for Astrophysics, Harvard University, 60 Garden Street, 

Cambridge, MA 02138 



> 

o 



and 



43 

6 

^ ■ Hsing-Wen Lin 3 



Institute of Astronomy, National Central University, 300 Jhunda Road, Jhungli, Taiwan 

32001 



\^Q' veres@ifa.hawaii.edu 
OV 

o; 

CN ■ Received ; accepted 



- 2 - 



ABSTRACT 



Asteroid detections in astronomical images may appear as trails due to a 
combination of their apparent rate of motion and exposure duration. Nearby 
asteroids in particular typically have high apparent rates of motion and accel- 
eration. Their recovery, especially on their discovery apparition, depends upon 
obtaining good astrometry from the trailed detections. We present an analytic 
function describing a trailed detection under the assumption of a Gaussian point 
spread function (PSF) and constant rate of motion. We have fit the function 
to both synthetic and real trailed asteroid detections from the Pan-STARRSl 
survey telescope to obtain accurate astrometry and photometry. For short trails 
our trailing function yields the same astrometric and photometry accuracy as a 
functionally simpler 2-d Gaussian but the latter underestimates the length of the 
trail — a parameter that can be important for measuring the object's rate of 
motion and assessing its cometary activity. For trails longer than about 10 pixels 
(> 3xPSF) our trail fitting provides ~ 3x better astrometric accuracy and up 
to 2 magnitudes improvement in the photometry. The trail fitting algorithm can 
be implemented at the source detection level for all detections to provide trail 
length and position angle that can be used to reduce the false tracklet rate. 

Subject headings: Solar system: Near Earth Asteroids, Data Analysis and Techniques 
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Introduction 



Since the advent of photographic plates asteroids have left their distinctive trails on 
astronomical images. These 'vermin of the sky' leave trailed detections because of the 
long time durations of the astronomical exposures and/or the asteroid's fast apparent 
rate of motion. An asteroid's trailing can be eliminated using non-sidereal tracking at the 
object's apparent rate of motion but this process simply transfers the trailing to the field 
stars while preserving the telescope system's point spread function (PSF) shape for the 
asteroid detection. Either way, the trailing spreads the total flux from the object over a 
larger area than the PSF, causing a reduction in the per unit area apparent magnitude 



and signal-to-noise, th ereby inducing a drop in the system's 



moving asteroids (e.g., 



Krugly 



2004 



Jedicke and Herron 



imiting magnitude for fast 



19971 ). Thus, dealing with trails 



on astronomical images is an unavoidable burden for astronomers who study the smaller 
objects in our solar system. 

In this work we introduce an analytic function describing a trail due to an object 
with an instantaneous Gaussian PSF moving at a constant rate of motion. We study the 
utility and efficiency of fitting the function to trailed detections using a large sample of 
both synthetic and real asteroids. A trail's morphological parameters as provided by the 
trail fitting can be used to significantly reduce the false tracklet detection rate which means 
that its use can deliver more discoveries of fast moving near-Earth objects (NEO), increase 
the upper rate of motion limit in automated surveys, and provide improved astrometry and 
photometry. 

Methods for identifying and characterizing trailed detections were introduced with 
the first CCD asteroid su rvey that used th e second order moments of a 2-d Gaussian to 



fit the trail (Spacewatch, 



Rabinowitz 



199lh . Lately, the characterization of trail length, 



astrometry and photometry has been used for artificial satellites and space debris studies 
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(e.g. iKouprianovl 120081 ; lLaas-Bourez et al.ll2009l ). This work deals specifically with asteroid 
trail fitting and characterization. 



Although more rea listic Lorentz and Moffat PSF profiles (jKouprianov 



Laas-Bourez et al. 



20081 : 



20091 ) have been used for the trail cores the analytical form of the 
trail equation has never been published. It is arguable whether a more complicated core 
is necessary given the prevalence of PSF changes during the exposure time combined 
with guiding and tracking irregularities. To reduce the computational load some earlier 
algorithms used a fixed or constrained trail length, orientation and/or brightness when 
fitting trails to known objects. Nowadays, modern all-sky surveys require robust image 
processing capable of reliable identification and characterization of transient detections 



including asteroid trails. 



The Panora mic Survey 



(Pan-STARRSl; 



Kaiser et al 



elesco pe and Rapid Response System's prototype telescope 



20101 )). the first of the next generation all-sky surveys, has 
submitted almost 3 million asteroid detections to the Minor Planet Center. Virtually 
all the reported astrometry and photometry are the result of a 2-d Gaussian fit to the 
detections using predicted (and fixed) PSF widths based on extrapolations from nearby 
field stars. While the resulting positions and flux are quite good for slow moving asteroids 
with PSF-like detection profiles we will show that the results can be improved dramatically 
with a true fit to the trails of fast moving asteroids like the NEOs. 

The trailing loss problem is illustrated in Fig. H] which shows how the error in the 
measured or reported magnitude drops as a function of the trail aspect (the ratio of the 
trail's length to its width). The figure clearly shows that the trailing loss can induce 
errors in the flux measurement that are much larger than the statistical uncertainty even 
for 'stubby' trails - those that are trailed only a few times more than their widths. A 
2-d symmetric Gaussian fit to the trail is clearly bad. In theory, an aperture magnitude 
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should work very well as long as the aperture is rectangular with its long axis aligned 
with the trail, of sufficient length to cover the entire trail, and with a width several times 
the PSF width. The catch-22 is that defining the aperture correctly implies knowing the 
correct dimension and orientation of the trail in the first place. The aperture magnitudes 
in Fig. H] made a simple assumption that the aperture has a dimension equal to 3 x the 
PSF width so it yields good results to trail aspects of ~3 but degrades quickly thereafter. 
Assuming larger apertures increases the likelihood of including background flux and thereby 
increasing the error on the derived flux. The real world Pan-STA RRSl telescope's Im age 



Processing Pipeline yields excellent photometry on point sources ( jSchlafly et al. 



20121 ) but 



fails on trailed sources because it does not implement a trailing fit. Presumably all the 
contemporary asteroid surveys suffer from this trailing loss at some level. 

Due to the worldwide interest in the NEOs, and their relatively high rate of motion 
that makes them difficult to recover if too long a time passes, NEO candidates are submitted 
as soon as possible on a per night basis as a set of > 2 detections known as a 'tracklet' 



(IKubica et al. 



20071 ). Subsequent inter-night linking of tracklets of the same object leads to 



derivation of the orbit and establishes whether the candidate is in fact a NEO. 

Tracklets are constructed by identifying a set of detections that are consistent with 
being the same object that have moved between exposures. The false tracklet rate is 
sensitive to the number of detections in the tracklet, the false detection rate, the maximum 
allowed rate of motion and acceleration, the astrometric and photometric error on each 
detection, the availability of the detections' morphological parameters, etc. In practice, the 
false tracklet rate is surprisingly high, usually due to systematic (rather than statistical) 



false detections. To reduce the false tra ck 
Processing System (MOPS; iDenneau et al. 



et ra te the Pan-STARRSl Moving Object 



20071 ) implements an upper limit on an object's 



rate of motion of 4 deg/day for tracklets containing > 3 detections and only 0.6 deg/day for 
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tracklets containing just 2 detections. Pan-STARRSl and contemporary surveys often rely 
on human observers to identify and verify trails and Pan-STARRSl is effectively limited 
to detecting objects with trails < 50 pixels due to the difficulty in identifying trails and 
the false tracklet explosion when attempting to link detections of objects that are moving 
quickly. Fitting all the detections to a trail fitting function helps in every way to reduce the 
false tracklet rate. 



2. The trailing functions 

2.1. PSF-convolution trail function 

An object moving with a constant apparent angular rate of motion u> (arcsec/second) 
in a CCD image with a detector pixel scale p (arcsec/pixel) in an exposure of time T 
(seconds) leaves a trail of length L = ooT/p pixels. When L is on the scale of the PSF 
width the source detection becomes non-PSF like — trailed. We approximate the trail as 
the convolution of an axisymmetric Gaussian PSF of width a moving at a constant rate in 
a direction x' that is rotated with respect to the image reference frame (x,y). The flux of 
the trail at any point in the (x',y') orthogonal coordinate system is then 

$ i f+ L / 2 

-L/2 



<f> 1 r+L/2 

f T (x',y') = b(x',y') + -^== / exp 

J-i y liter 1 J -l/2 



1 {{x'-iY + wY} 



di (i) 



la 1 

where $ is the total photometric flux in the trail and b(x',y') is the background flux at the 
same point. The flux in any pixel is strictly the integral of g(x',y') over the bounds 
of the pixel but in images that are not under-sampled the function changes slowly enough 
across the pixel that we can assume the pixel flux is simply the flux at the center of the 
pixel. 

This equation can be rewritten using the error function erf (z) = ^= e~* 2 dt, rotating 
the trail to the image reference frame through an angle 9 with respect to the +x-axis, and 
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translating the trail's centroid to (xo,yo) such that 



x' = (x — xq) cos 9 — (y — ?/o ) sin 6* 
y' — (x — Xq) sin # + (7/ — t/ ) cos & 



(2) 



yielding the trail equation 



$ 1 

/tO, y) = b{x, y) + — - — f== exp 



^2aV27r 



{{x - x ) sinfl + (7/ - y ) cos9)f 
2a 2 



erf 



(x — Xq) cos 9 + (y — t/o ) s i n ^ + 



(7 



v/2 



erf 



(x — x ) cos 9 + (y — 7/o ) srn ^ — -^/2 



V2 



(3) 



It is well known that an axisymmetric Gaussian profile underestimates the flux in real 
PSFs and other more elaborate PSF models such as the Lorentz or Moffat functional forms 
could be used (or e ven empirical functions based on the realized PSFs within an image (e.g. 



Trujillo et al. 



20011 ) ). Any core PSF function can be convolved with a line to provide a trail 
function though most will of necessity need to be solved numerically. However, the utility of 
doing so is questionable given that in real images the PSF shape and flux are typically not 
constant in time — the object leaving a trail may itself be changing its intrinsic flux; the 
seeing can change during the exposure as may the optical system's transmission function 
due to system flexure; wind buffeting, tracking and guiding can change the apparent path 
of the object on the CCD or its rate of motion across the pixels; etc. 



2.2. 2-d Gaussian trail function 



Elongated detections have classically been recognized when the major axis 2 nd moment 
significantly exceeds the 2 nd minor axis moment in a fit to a 2-d Gaussian: 



faW, y') = K x ' v') + t, — — ex P 



X 



r> 



y 



2<£ K 



(4) 
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with widths (ov, a y r). Rotating the function using eq. [2] yields 



fa{x,y) = b(x,y) + 



2n<j x (jy 



exp 



(x-x ) 2 (y-yof 



2a 2 



2a 2 , 

y 



(5) 



Fitting a trail to this functional form provides estimates for the centroid position, second 
order moments, background, and total flux — but not the trail length. 

For the 2-d Gaussian we use the major axis 2 nd moment as a proxy for the trail length: 



L ~ L g — a x , 




2x 



f 1/2L $ 1 



1/2L 



L2na 2 



exp 



2a 2 



\ {(x -x ) 2 + y' 



dxo dx dy . (6) 



Integrating equation ([6]) and dropping the prime notation we find that: 



L g = Jl2 (a 2 -a 2 ) 



(7) 



2.3. Acceleration &; phase angle effects 

Ignoring the topocentric acceleration (a) in the trail fit causes an astrometric error at 
the mid-time of the exposure, the typical time at which astrometry is reported, of a T 2 /8 
if all the acceleration is parallel to the direction of motion (along-track). Taking as an 
upper limit a good astrometric error on reported observations of 0.1" than a T 2 /8 < 0.1" 
or a < 1"/ T 2 to guarantee that the acceleration does not induce significant error in the fit. 
For point of reference, typical exposure times for contemporary asteroid surveys are on the 
order of 100 s so that accelerations > 10~ 4 arcsec/s 2 induce significant astrometric error. 

The fastest topocentric accelerations are observed for objects that are very close to 
Earth so we examined the apparent rate of motion and acceleration of known asteroids at 
the moment of their closest approach within 10 lunar distances between 1900 A.D. and 
2200 A.D. (see figs. H]-|3]). Of the 601 close approaches only about 2% have accelerations 
on the order of the canonical limit of 10~ 4 arcsec/s 2 . Objects on heliocentric orbits exhibit 
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these accelerations when they are within about one lunar distance (LD), more typically 
about 0.5 LDs. We note that these objects are moving so quickly across the image plane 
that they leave very long trails at typical pixel scales — since there is no advantage to long 
exposures for these trails (the flux per unit trail length is constant) the solution is simply 
to use shorter exposures to eliminate the astrometric impact of neglecting the objects' 
acceleration. 

Another possible issue for the astrometric fit to a trail is the effect of changing 
brightness of the object during the coarse of a single exposure due to the changing phase 
angle. Figure [2] shows that the fastest moving objects at closest approach to Earth 
change their brightness at >0.1 mmag/s so a 100 s exposure would experience >0.01 mag 
change in brightness with a mean of ~0.03 mags. This effect is much smaller than typical 
photometric uncertainty on even slow moving asteroids and comparable to photometric 
errors induced by variable seeing conditions and sky catalogs. The small asteroids that 
are typically observed within 1 LD are also rotating on time scales comparable to typical 
exposure times and the rotation a mplitudes are much larger than the phase angle effects 



( iHergenrother and Whiteley 



20111 ) . Furthermore, linear changes in the brightness will 
probably not induce large errors in the fitted centroid — only non-linear changes in the 
flux that occur at an even smaller level. We conclude that the astrometric impact of 
changing phase angle is extremely small and, in the worst cases as described above, can be 
ameliorated simply by using shorter exposure times. 

Thus, while this work is strictly applicable only to zero acceleration trails with no 
change in brightness during the exposure, the vast majority of objects exhibit negligible 
acceleration and brightness change within the exposure time and our technique can be 
applied without astrometric impact. The method is applicable as long as a T 2 /8 is less 
than the typical astrometric uncertainty which can almost always be achieved simply by 
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decreasing the exposure time. 



3. Trail fitting performance 

We measured the performance of the trail fitting functions on synthetic and real trails in 
Pan-STARRSl images. The synthetic trails were created to mimic the real Pan-STARRSl 
trails but allowed us to control the exact flux, position, length and orientation of the trails 
to determine how well the fitting procedure reproduce the generated values. 

We fit the trails to the functions using the IDL procedure mpfit2dfun0 that employes 
the Levenberg-Marquar dt least-squares fitting technique to minimize the variance between 



the trail and the model (ILevenberg 



19441 : 



Marquardt 



1963|). 



We imposed several constraints on the fitting procedure to reduce the computation 
time and reduce the likelihood of non-sensical results. We set the maximum number of 
iterations to 50 after determining that all good fits were identified by the algorithm in < 50 
iterations and imposed limits on the ranges of the fit parameters — a, L, xq, y^. The trail 
width a is constrained to be within the range of stellar PSF widths from all Pan-STARRSl 
images (1.5 pix < o <4.5 pix), L must be smaller than the image dimension, and the 
centroid position (xq, yo) must be within 20 pixels of the center of the image. The last 
constraint was implemented because the image 'source' detected by the source detection 
algorithm was placed at the image center. 

The resulting fit returned the best values of the 7 fit parameters, x = (b, $, L, a, 9, xq, yo), 
and their associated uncertainties, 5x = (5b, £$, SL, 8a, 59, 5x , 5y ). We then calculated the 



1 Markwardt IDL library, http://www.physics.wisc.edu/~craigm/idl 
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X 



2 of the fit in a 'trail aperture': 



aperture 



[ffit(i,r, x) - f(i,j)] 2 
8fw(i,j; x,5x) 



* 2 = £ 



8 



where f(i,j) is the actual flux in the pixel (i, j), ff it (i,j) is the value of the fit function at 
pixel and Sff it (i,j) is the propagated uncertainty in the fit at the same pixel. The 
'trail aperture' over which the summation is performed has a rectangular shape with center 
at the fitted centroid position (x ,y ), length of L + 6a, width of 3a, and the long axis 
of the rectangular aperture is aligned with the major axis of the fitted trail. The reduced 
chi-squared is xled = X 2 /(N — 7 — 1) where the denominator is the number of degrees of 
freedom in the fit when N is the number of the pixels in the aperture. We remind the 
reader that Xred ~ 1 f° r a g°°d fit while much larger or smaller values indicate a bad fit. 

In the following sections we quantify the trail fitting algorithms' performance as a 
function of the trails' statistical significance: S/N = ^g +B - The signal S is the integrated 
flux $ of the source as measured by the fit or as generated for synthetic trails. The 
background B is the integrated total background flux in the trail aperture including readout 
noise, dark current, sky, and other sources of flux that are not due to the trailed object. 



Synthetic trails were created using eq. |3j In most astronomical images the background 
is dominated by the sky and over a small region near a trail the sky is generally flat after 
correcting for any flat-field issues. Thus, in the absence of noise the flux at a given pixel is 
fsyn{hj) = f{hj] x) + B where B is a constant background level. In the presence of random 



ran(l) is a random number from the normal distribution with centroid and width 1. 

Figure [5] provides a few examples of synthetic trails where the only difference between 



3.1. Generating synthetic trails 



statistical noise the signal at the pixel is f' syn (i,j) = f S yn{ 



+ ran(l)y7 S2/n (i, j) where 
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them is their S/N. The faintest trail at S/N = 7 is only visible as a faint smudge and 
the eye is guided to it because we know that the trail lies at the center of that sub-image. 
Detecting these trails with automated software is not simple — measuring their properties 
once detected is the purpose of this work. 

Figure E] illustrates the difference between the 'true' trail's profile (eq. [3]) and a fit to 
a 2-d Gaussian (eq. [6]). It is clear that in this example the latter can provide reasonable 
estimates of the trail's relevant parameters — its flux, centroid and length (perhaps after 
correction with a simple multiplicative length-based factor). We will show later that the 2-d 
Gaussian's utility breaks down for longer trails but its qualitative match to the generated 
synthetic trails degrades even sooner as shown in figure [7J 



3.2. Synthetic trail fitting 

We measured the trail fitting performance on 30,000 synthetic trails generated with 
a flat distribution in length (10 < L < 50), orientation (-90° < 9 < +90°), width 
(1.5 pixels < a < 4.0 pixels), and 3 < S/N < 65. The range in trail length corresponds to 
apparent motions of 1.3 to 13.3 deg/day in the Pan-STARRSl Solar System Survey mode 
with 45 sec exposure times and a pixel scale of 0.25 arcsec/pixel. The centroids of the 
synthetic trails were fixed at the center of the small synthetic images. 

Figure [S] shows that, as expected, both astrometric error (A x .) and uncertainty (a x ) 
become smaller as S/N increases and a x > A x . At marginal S/N the error is on the order 
of 2 pixels. With this metric and for these synthetic trails the 2-d Gaussian fit yields 
astrometric errors that range from about ~ 1.5x to ~ 2.8x higher than the trail fitted 
errors for 10 < S/N < 50. While the astrometric errors are much improved with the trail 
fitting the reported uncertainties are essentially identical. Thus, the 2-d Gaussian fitting 
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misrepresents the actual error. 

The utility of the trail fitting function is evident in fig. M which illustrates how the 
astrometric error and uncertainty evolve with the trail 'aspect': L/a, i.e. the length of the 
trail relative to the PSF. A detection with L/a ~ 1 appears untrailed while detections with 
L/a > 5 are clearly trailed to eye. The figure shows that as a trail's aspect increases both 
the astrometric error and uncertainty increase until they plateau for L/a > 10. The plateau 
values are ~50-100% higher for the 2-d Gaussian fit compared to the trail fitting function. 

A problem with trailed asteroid detections is the 'loss' of S/N as the object becomes 
trailed. There are two ways to think of this effect: 1) the peak signal per pixel decreases 
as 1/L while the per pixel noise remains roughly constant or 2) the total signal in the trail 
remains constant as the noise increases because there are more pixels 'under' the trail. 
The former interpretation is a trailing loss because it affects the ability to detect the trail 
using peak-pixel detection algorithms. The latter interpretation affects the ability of more 
sophisticated algorithms that attempt to identify trails by integrating the flux along lines. 
In each case the ability to detect a trailed asteroid decreases due to trailing loss. The 
amount of trailing loss is not the same in the two classes of algorithms and the implications 
are too complex to consider in this work that only discusses fitting the trails once they have 
been identified. 

Once the trail has been detected it is important to correctly measure the trail's flux 
and fig. [10] illustrates the results of the flux determination using the trail and 2-d Gaussian 
fitting methods. It is clear that the trail fitting is superior to the 2-d Gaussian fit but the 
figure under-represents the improvements with the trail fitting because it includes trails of 
all lengths at each total S/N. At small total S/N for longer trails the flux is spread over 
many pixels so that the per pixel S/N is much smaller (the 'trailing losses' are described 
in more detail later in this section). In this situation, by far the most common since the 
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number of trails increases exponentially with decreasing S/N, the trail fitting algorithm 
dramatically outperforms the 2-d Gaussian method. 

Figure [11] provides a comparison to the aperture measured photometry using apertures 
specific to the trail and fit type. i.e. the trail fit uses a rectangular aperture as described in 
§121 and the 2-d Gaussian aperture uses a similar rectangular apertue but with L — > L g and 
u — y Ug. We see that trail photometry using apertures based on the generated trail lengths 
and widths reproduces the expected flux with zero error as expected. An aperture based 
on the trail- fit length and width performs very well for S/N > 5 while a flux measurement 
based on the 2-d Gaussian fit aperture only has approximately zero error for S/N > 25. 
Once again, since most trails have small S/N- weighted error metric would clearly show that 
the trail-fitting technique offers a major advantage over the 2-d Gaussian fit. 

While the primary purpose of this work was to illustrate the astrometric and 
photometric improvements when using our PSF-convolution trail fitting the remaining 



two parameters, the trail length and orientation, are i mportant in 



inking individual 



20071 ). This is a 



detections on the same night together into 'tracklets' (IKubica et al. 
common combinatoric explosion when asteroid surveys attempt to link detections of fast 
moving asteroids — the search area increases as the square of the upper limit on the rate 
of motion thereby squaring the number of candidate matching detections. One might think 
that the astrometric and photometric requirements would be sufficient to eliminate false 
linkages between the detections but, surprisingly, there are often far more false detections 
than statistics would predict leading to a disturbing number of false tracklets. The typical 
solution is to decrease the upper rate limit to keep the false tracklet rate manageable and 
to increase the number of required detections for a valid tracklet. 

Another solution to the tracklet formation problem is to use more information - 
larger areas must be searched as the upper limit on the rate of motion is increased but the 
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faster moving objects will also leave trails rather than point source detections. Those trails 
will typically be of approximately the same length (assuming constant exposure times) and 
will have roughly the same orientations. Thus, a good trail fitting algorithm on candidate 
detections can dramatically reduce the false tracklet formation rate. 

The ability to correctly measure the length of a trail is illustrated in fig. [T^] in 
terms of the derived trail aspect error and uncertainty as a function of S/N. The length 
derived with the 2-d Gaussian is always underestimated by > 10% for S/N > 5 while the 
PSF- convolution trail fitting provides an error of < 3% in the same S/N range. This is not 
unsurprising as, by its very nature, the 2-d Gaussian should always underestimate the trail 
length. The trail aspect uncertainty returned by the 2-d Guassian is always underestimated 
by about 10% while the trail fit uncertainty is always slightly higher than the error as 
desired. 

Similarly, fig. [13] demonstrates that the error in orientation does not depend on the 
S/N but that the PSF-convolution method consistently performs better than the 2-d 
Guassian fits. Since the measured trail lengths/aspects and orientations can be used to 
reduce the false tracklet rate it is important to use a functional form that returns more 
accurate values and the PSF-convolution trail fit is clearly superior. 

4. Real trail fitting 

We have shown above that our PSF-convolution trail function performs well on fitting 
synthetic trails but the true litmus test is how it performs on real detections in real images 
with all their systematic noise problems, dead and bright pixels, edges and gaps, saturated 
stars and diffraction spikes, ghost images, cosmic rays, etc. 

To test our algorithm on real trails we employed asteroids detected by the 
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Pan-STARRSl telescope and identified by MOPS. The MOPS database contains hundreds 
of thousands of known and unknown asteroids of which several thousand are trailed. 
MOPS stores 200x200 pixel 'postage stamps' of all Pan-STARRSl detections that were 
incorporated into tracklets. The postage stamp FITS file is the differen ce of two images 



20071 ) and is therefore 



of the same portion of the sky exposed at different times (jLupton 
already background subtracted. 

Figure [H] provides examples of a couple real postage stamps — one ideal and one 
typical. In the ideal example our trail fitting worked perfectly but even for the more typical 
case the PSF-convolution trail fit provides a good reconstruction of the visible part of the 
detection. Even though the trail has been truncated there is still useful information in its 
orientation and length to assist in linkning similar detections on a night into tracklets. In 
cases where the trail is truncated our trail fitting algorithm raises a flag that there are 
masked pixels in the aperture. In these cases the astrometry and photometry are also 
flagged as suspect. 

To measure the astrometric error for real asteroid trails we selected the 1,000 longest 
trails of numbered asteroid detections from the MOPS database that were submitted to the 
MPC. These trails have an average length of 15 pixels and the longest automatically detected 
trail in the sample is 25 pixels long. e.g. an average trail aspect of ~3. For comparison, 
the average length of the 1,000 longest submitted tracklets irrespective of whether they are 
known is 18 pixels with the longest being 43 pixels. Since the Pan-STARRSl survey has not 
emphasized the detection of trailed, fast moving objects, even longer trails exist in the data 
that were not detected by the IPP's source detection algorithms. For instance, we manually 
identified 114 pixel long trailed detections of 2012 LZi during its close approach to the 
Earth. These and other trails will be identified and measured as the data are reprocessed 
but the objects will be impossible to recover. On the other hand, their astrometry and 
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photometry will be useful for future precoveries of the objects by Pan-STARRSl and other 
surveys. 

We calculated the expected position for each of the 1,000 trailed detections of 



20091) that 



numbered asteroids at the epoch of observation using OpenOrb (IGranvik et al. 
uses DE406 as input for the perturbing planets. Since the numbered asteroids have accurate 
orbits we can use their calculated positions to measure the astrometric error of our trail 
fitting method. Figure [T5l shows that the PSF-convolution trail fitting provides a 2-3 x 
improvement in the realized astrometric error over a 2-d Gaussian fit for these trails. 

It is clear that figure [15] underestimates the improvement in astrometric error in the 
trail fitting because we have tested the algorithm only on trails for numbered asteroids. 
These objects are on average larger, more distant, and leave shorter trails (smaller trail 
aspects) on the image than the smaller, closer, unnumbered and newly discovered asteroid 
trails. Reducing the astrometric (and photometric) error on newly discovered objects 
will have a dramatic impact on the ability to recover the objects shortly after discovery 
by 1) decreasing the sky-plane ephemeris uncertainty and 2) providing better magnitude 
estimates to improve target selection by the followup observatories, i.e. to better match the 
sites' detection capabilities to the target's brightness. 

As discussed above, false linkages can be dramatically reduced by requiring that the 
detections within a tracklet have the same orientation (position angle; PA) and roughly the 
same lengths. The fitted length of real trails is not as robust a filter as the PA because real 
trails often intersect chip gaps, boundaries, diffraction spikes, etc., that tend to truncate the 
measured trail length. In practise, the PA filter constraint must be trail length dependent. 

Figures [TBI and [TTI show the relative error in both measures in our 1,000 trail sample 
and generally mimic the results from our studies with synthetic trails as shown in figs. [12] 
and [13] The PSF-convolution trail fit provides lengths (or equivalently, aspects) consistent 
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with the ephemeris prediction but the 2-d Gaussian derived lengths are underestimated by 
< 20%. The PSF-convolution trail fit derived orientation matches the predicted position 
angle for all aspect ratios but the 2-d Gaussian-derived orientations show a strange and 
unexplained error of up to 10% for moderate trail aspect ratios. 

One of the major problems with trail identification and characterization is computing 
time. Even though the human eye and mind are tuned to identify these features writing 
a computer program to do so in astronomical images and at small S/N has proven to 
be challenging. Fitting our trail function to a source detection that has already been 
identified is less time consuming but still ~ 3x more computationally expensive than simple 
PSF fitting. For example, Pan-STARRSl identifies so many false detections that it is 
computationally impossible with the current hardware resources to fit all identified source 
detections to the PSF-convolution trail function even though doing so would dramatically 
reduce the false detection rate. 



Instead, MOPS currentl y links detections 



motion between detections ( IMilani et al. 



Dased purely on the astrometry and consistent 



20121 ) and only fits the PSF-convolution trailing 



function to those detections vetted and selected by a human observer. Along with the 
fitted trail's \ 2 we find that we can select for good trails simply by limiting the number 
of iterations in the fit (see fig JTBl . The fits to fully 99% of real asteroid detections (both 
trailed and PSF-like) converged in under 25 iterations while 69% of false detections do not 
converge. (The number of iterations needed for convergence was tested on a sample of 1,000 
trailed, 1,000 PSF-like and 1,000 false detections from Pan-STARRSl.) This is a powerful 
computational savings since the number of false detections far outnumbers the number of 
real detections. 
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4.1. False detection reduction 

Trail fitting is effective at eliminating false detections and even better at eliminating 
false tracklets. In Pan-STARRSl images one of the major problems is false detections 
along diffraction spikes as shown in fig. [TJO Statistical fluctuations along the length of the 
diffraction spikes trigger multiple instances of the PSF-detection algorithm. These false 
detections can be eliminated by subsequently fitting them to our trailing function because 
the trail fit returns an unreasonably large RMS and x 2 an d position angles that are aligned 
with the expected PA for diffraction spikes in the image. 

We measured the trail fitting's false detection reduction capability using postage 
stamps on ~34,000 detections identified by Pan-STARRSl's PSF-detection algorithm. The 
maximum number of iterations was set to 25 (as per figure [TBI and we required that the 
fit's reduced x 2 < 3. These cuts kept ~96% of real detections and rejected 98% of the false 
detections. 

Tracklet rejection was based on cuts suggested by figures [16] and [T7] e.g. detections 
within a tracklet must have lengths and orientations consistent to within 10% and, as 
discussed above, the fits can not be compromised by contacting masked detector areas. We 
achieved a false tracklet rejection rate of 99% while the fraction of rejected real tracklets 
was only ~1%. The tracklet rejection rate can be higher than the detection rejection rate 
because the same false detections can appear in multiple tracklets — each of which is 
rejected by identifying the errant detection. 

5. Conclusions 

We provide the analytic form for an asteroid trailing function assuming that a 
symmetrical Gaussian PSF is moving at constant speed during an exposure. We then 
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demonstrate that this function provides accurate photometry and astrometry — much 
better than the values derived using 2-d Gaussian fits to the trails and much better than 
the values currently being provided by the Pan-STARRSl survey. 

The trailing function provides astrometry with > 3x smaller errors than are typically 
reported by the surveys. The photometry with the trail fitting improves by an even larger 
factor. It is clear that the trail fitting function should be implemented whenever possible 
- even when fitting stationary objects since they are nothing but trails of length L = 0. 
We note that there are extra benefits to using the trailing function for all image sources 
because the length and orientation can be used to diagnose optical and focus effects within 
the image plane. 

The trail fitting morphological parameters can eliminate on the order of 99% of 
false tracklets by requiring that detections within the tracklet have the same lengths, 
orientations, and fluxes. The derived trail lengths and position angles can also help to 
identify false detections that are saturated streaks and diffraction spikes. 

The trail fitting computation time might be an issue because it can be significantly 
longer than for the 2-d Gaussian. A hybrid method that will be implemented by 
Pan-STARRSl is to apply the trail fitting only to detections with large second order 
moments and/or excessive flux in the aperture that might indicate a trail. 
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Fig. 1. — Geocentric distance and apparent rate of motion of known asteroids at the moment 
of closest Earth approach between the years 1900 and 2200 A.D. 
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Fig. 2. — Apparent rate of motion, acceleration, and rate of change in apparent magnitude 
of known asteroids at the moment of closest Earth approach between the years 1900 and 
2200 A.D. (Left) Apparent rate of motion vs. apparent acceleration (Right) Rate of change 
in apparent magnitude vs. apparent rate of motion. 
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Fig. 3. — Along- and cross-trail components of the acceleration of known asteroids at the 
moment of closest Earth approach between the years 1900 and 2200 A.D. Rates of accelera- 
tion above and to the right of the dashed curves are large enough to induce 0.1" astrometric 
error in the fitted trail position at the specified exposure times of 30, 60 and 100 seconds. 
i.e. the technique described in this paper does not induce significant astrometric error for 
rates of acceleration in the lower left corner of the figure. 
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Fig. 4. — Trailing 'loss' as a function of the trail's aspect ratio (the ratio of its length to 
width). A perfect algorithm for measuring the magnitude of a trailed detection would provide 
zero magnitude loss, (solid) The average derived loss in magnitude for synthetic trails fit to 
2-d symmetric Gaussians with width equal to the trail width i.e. equal to the PSF. (dashed) 
The average derived loss in magnitude for synthetic trails using square apertures with sides 
equal to 3x the trail width, (data points) The difference in magnitude between the reported 
Pan-STARRSl instrumental magnitudes and the correct value for the 1,000 longest real 
trails for known asteroids. 
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Fig. 5. — Three synthetic images of synthetic trails with S/N— 7, 25, 85. All the trails have 
L = 30 pixels, cr=2 pixels and 6 = n/5 = 36 deg. 
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Fig. 6. — Profile of a synthetic trail (L = 30, a = 2, S/N=70) including noise (solid) along its 
major axis along with the generated trail's shape (dashed) and the result of a 2-d Gaussian 
fit (dotted). 
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Fig. 7. — (left) Generated synthetic trail (L=40, a = 2.5, 9 = vr/4), the same trail recon- 
structed (center) by the trail fit and the same trail reconstructed with a 2-d Gaussian fit 
(right). 




Fig. 8. — (Left) Moving average of the astrometric error for the PSF-convolution and 2-d 
Gaussian trail fits as a function of the trail's S/N. The error bars are the standard error on 
the mean. (Right) Same as on the left but for the reported uncertainty. 
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Fig. 9. — (Left) Moving average of the astrometric error for the PSF-convolution and 2-d 
Gaussian trail fits as a function of trail aspect (L/a) for synthetic trails with 3 < S/N < 65. 
The error bars are the standard error on the mean. (Right) Same as on the left but for the 
reported uncertainty. 
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Fig. 10. — (Left) Fractional error and (right) uncertainty of the flux derived with the trail 
and 2-d Gaussian fitting as a function of the trail's total S/N. The results include trails of 
all lengths at the same S/N as described in the text. 
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Fig. 11. — Fractional error of the measured trail flux in rectangular apertures with size 
and orientation determined by the two types of trail fitting and using the generated (input) 
parameters versus total S/N. The results include trails of all lengths at the same S/N. 
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Fig. 12. — (Left) Moving average of the fractional error on the trail aspect (L/cr) for the PSF- 
convolution and 2-d Gaussian trail fits as a function of S/N. The error bars are the standard 
error on the mean. The trail aspect for these synthetic trails ranges from 3 < L/a < 35. 
(Right) Same as on the left but for the reported uncertainty. 
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Fig. 13. — (Left) Moving average of the fractional error on the trail orientation for the PSF- 
convolution and 2-d Gaussian trail fits as a function of S/N. The error bars are the standard 
error on the mean. The trail aspect for these synthetic trails ranges from 3 < L/a < 35. 
(Right) Same as on the left but for the reported uncertainty. 



Fig. 14. — (left) Real Pan-STARRSl difference images of trailed asteroids and (right) their 
reconstructed trails from our PSF-convolution trail fitting, (top) An almost ideal case with 
a nicely trailed detection in the center and no obvious image artifacts and (bottom) a more 
typical case where the trailed detection passes through the masked pixels of a badly sub- 
tracted bright star with diffraction spikes, 'streaks', and a gap between CCD chips. Masked 
or non-existent pixels are in grey. 
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Fig. 15. — (solid) Astrometric error after trail fitting for the 1,000 longest trailed detections 
of numbered asteroids that were submitted to the MPC by Pan-STARRSland (dashed) for 
the submitted astrometry from the Pan-STARRSl's PSF fitting algorithm. For compari- 
son, the average P an-STARRSl astrometric error for untrailed detections is ~ 0.13 arcsec 



(IMilani et al. 



20121 ). 
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Fig. 16. — Fractional error in the trail aspect L/a for the 1,000 longest trails for known 
asteroids for both the (solid) trail fitting and (dashed) 2-d Gaussian fits. 
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Fig. 17. — Fractional error in the trail orientation 9 for the 1,000 longest trails for known 
asteroids for both the (solid) trail fitting and (dashed) 2-d Gaussian fits. 
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Fig. 18. — Number of iterations needed for fit convergence for trailed, PSF-like and false 
detections. 
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Fig. 19. — Four examples of (top) diffraction spikes in difference images and (bottom) their 
corresponding trail fits. The difference images show both the positive and negative diffraction 
spikes from the same stars but with field rotation between the two images. The trails are 
only fit to the positive diffraction spikes. The thick gray bands in the difference images 
represent chip gaps. 



